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ABSTRACT 

The knowledge of accurate stellar parameters is paramount in several fields of 
stellar astrophysics, particularly in the study of extrasolar planets, where often the 
star is the only visible component and therefore used to infer the planet's fundamental 
parameters. Another important aspect of the analysis of planetary systems is the stel- 
lar activity and the possible star-planet interaction. Here we present a self-consistent 
abundance analysis of the planet hosting star WASP- 12 and a high-precision search 
for a structured stellar magnetic field on the basis of spectropolarimetric observations 
obtained with the ESPaDOnS spectropolarimeter. Our results show that the star does 
not have a structured magnetic field, and that the obtained fundamental parameters 
are in good agreement with what was previously published. In addition we derive 
improved constraints on the stellar age (1.0-2.65 Gyr), mass (1.23-1.49 M/Mq), and 
distance (295-465 pc). WASP- 12 is an ideal object to look for pollution signatures 
in the stellar atmosphere. We analyse the WASP- 12 abundances as a function of the 
condensation temperature and compare them with those published by several other 
authors on planet hosting and non-planet hosting stars. We find hints of atmospheric 
pollution in WASP-12's photosphere, but are unable to reach firm conclusions with 
our present data. We conclude that a differential analysis based on WASP- 12 twins 
will probably clarify if an atmospheric pollution is present, the nature of this pollution 
and its implications in the planet formation and evolution. We attempt also the direct 
detection of the circumstellar disk through infrared excess, but without success. 

Subject headings: stars: individual (WASP-12) — stars: abundances — stars: magnetic field 
— stars: fundamental parameters 



1. Introduction 



One of the biggest surprises in the exoplanet field was the discovery of gas giant planets 
orbiting very close to their host star. These hot Jupiter planets represent one extreme of the 
Galaxy's population of planets, and they provide important constraints to guide our nascent ideas 
about the formation an d evolution of plane tary systems. The probability of transit for close-in 



giant planets is ~ 10% (|Seager et al 



2000|) . and through the analysis of transit light curves the 



ratio of the stellar and planetary radii can be deduced. Through the significant uncertainties in 
the mass and radius of any particular s tar, our characterisation of exoplanets is limited by that of 



their host stars (e.g. 



Southworth 



I2009|) . For this reason, it is important to directly measure the 
properties of planet hosting stars, particularly in cases where we expect the presence of planets 
may have influenced the properties of the star through star-planet interactions. 

One of the most extreme hot Jupiter e xoplanets is WAS P- 12 b, a gas giant planet orbiting 



only 0.023 AU from a late F-type host star (|Hebb et al 



2009). WASP-12b's orbit is, therefore. 



only about 1.5 stellar diameters from the photosphere of the star. At such proximity, interactions 
between the star and the planet must occur. Near-UV observations of WASP- 12 covering the 
wavelengths of many resonance lines re veal that WASP- 12 b is surrounded by an exosphere 



which appears to overfill its Roche lobe (IFossati et al. 



20101) . This exospheric gas may be the 



consequence o f tidal disruption of the planet's convective envelope as recently suggested by 



Li et al 



(120101), but it could also be entrained material from the stellar corona. 



A planet or biting as close as WASP- 12 b might be expec ted to interact magnetically with 



its host star, c.f. 'Shkolnik et al. 



(12003b, 



Shkolnik et al. 



(l2005h . and 



Shkolnik et al. 



(12008b. A 



first step to search for si ich interactions is to detect and quantify the stellar magnetic field by 



spectropolarimetry, e.g. 



Fares et al. 



(120101) . The presence of a magnetic field belonging to 
WASP- 12 would provide a precious piece of inforr nation needed t o establish which mechanism 
controls the structure and the evolution of the disk tLai et al.i.2010b . 
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The solar system's giant planets have enhanced metal abundances relative to the Sun (IGuillot 



2005|) and high atmospheric metal abundances have been suggested as a contributi ng factor in 



the infl ated radii of planets such as WASP-12b ([Burrows et al. 



20071). Intriguingly, 



Fossati et al 



(|2010() detected a wealth of metallic atoms and ions in the exosp heric gas 



Li et al. 



surrou nding WASP- 1 2 b. 



(l2010h . this could lead to 



If this gas is indeed accreting onto the host star as suggested by 
abundance anomalies in the photosphere of WASP- 12. Since WASP- 12 is expected to have a very 
shallow surface convection zone, any accreted gas will remain close to the surface rendering any 
pollution of the surface composition relatively easy to detect. 

In this paper, we report on spectropolarimetric observations of WASP- 12 which we use 
to probe the stellar magnetic field, fundamental parameters and abundance pattern of the 
star. In Sect. [2] we describe our observations and data reduction. Sect. [3] and Sect. |4] provide 
a description of the adopted model atmosphere including methods and results of the stellar 
parameter determination and abundance analysis. In Sect. [5] we provide the results of our stellar 
magnetic field search. Our results are finally discussed in Sect. [6l while in Sect. |7] we gather our 
conclusions. 

2. Observations and data reduction 



We observed WASP- 12 using the ESPaDOnS (Echelle SpectroPolarimetric Device for 
ObservatioNs of Stars) spectropolarimeter at the Canada-France-Hawaii Telescope (CFHT) on the 
3rd and 5th of January 2010. The observations were performed in "polarimetric" mode. 

ESPaDOnS consists of a table-top cross-dispersed echelle spectrograph fed via a double 
optical fiber directly from a Cassegrain-mounted polarisation analysis module. Beside the natural 
intensity /, in polarimetric mode the instrument can acquire a Stokes V (or Q or U) profile 
throughout the spectral range 3700-10400 A with a resolving power of about 65 000. A complete 



polari metric observation consists of a sequence of 4 sub-exposures (|Donati et al .11 19971 : 1 Wade et al 



2000 ) 



Each of the four sub-exposures was 1290 seconds long, with a total amount of integration 
time of 1.5hrs, eac h night. The spectra were reduced using the Libre-ESpRlT package^ 



(Donatiet al. 



20071). The Stokes / spectra have a signal-to-noise ratio (SNR) per pixel of ~ 126 



and ~158 in the continuum, on the first and second nights respectively; both values calculated 
at 5000 A. To increase the SNR of the Stokes / spectrum we averaged the two available spectra, 
obtaining a single spectrum with a SNR of ~200, normalised by fitting a low order polynomial to 
carefully selected continuum points. 

The eff'ective temperature (Tes) was determined (see Sect. S]) from our ESPaDOnS data and 
two spectra of WASP- 12 obtained at the Isaac Newton Telescope (INT) with the Intermediate 
Dispersion Spectrograph ( IDS). The spectra cover the region of the Hor line with a spectral 



resolution of i?=8 000 (see 



Hebb et al 



20091 for more details). 



W e also used Near-UV observations obtai ned with the HST Cosmic Origin Spectrograph 

Lo|;l 



(COS) (Gree n et all 



201( 



Osterman et al. 



20101) for the analysis of the spectral energy distribution. 



The spectra, calibrated in flux, cover three non-contiguous wavelength rang es in the Nea r -UV 



with a resolution of i? ~20000. These observations are described in detail in 



Fossati et al. 



( 20101) 



3. The model atmosphere 



To compute model atmospheres o f WASP- 12 we employed the LLmodels stellar model 



atmosphere code dShulyak et al 



2004) . For all the calculations Local Thermodynamical 



Equilibrium (LTE) and plane-parallel geometry were assumed. We used the VALD database 



'www . ast . obs-mip . £r/pro j ets/espadons/espadons . html 
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(IPiskunov et al. 



1995 



Kupka et al 



19991 : 



Rvabchikova et al 



1999 ) as a source of atomic line 



parameters for opacity calculations. The recent VALD compilation contains information for 
about 6.6 X 10' atomic transitions, most of them coming from the latest t heoretical calculations 



1991 



J 



performed b y R. Kuruc^tj. Convection was implernented according to the 



1992h model of convection (see 



Heiter et al 



Canuto & Mazzitelli 



2OO2L for more details). 



4. Fundamental parameters and abundance analysis 



Hebb et all (|2009h derived the fundamental parameters of WASP- 12 from the analysis of low 
and mid-resolution spectra, obtaining Teff =6300+ 150 K, log ^=4.38+0. 10, and adopting a value of 
0.85 km s'^ for the microturbulence velocity (i/mk)- We used these values as our starting point in 
an iterative process to gradually improve the parameters using different spectroscopic indicators. 
In our analysis, every time any of T^s, log^, fmic or abundances changed during the iteration 
process, we recalculated a new model with the implementation of the last measured quantities. 
Similarly the derived abundances were treated iteratively: while the results of the abundance 
analysis depend upon the assumed model atmosphere, the atmospheric temperature-pressure 
structure itself depends upon the adopted abundances, so we recalculated the model atmosphere 
every time abundances were changed, even if the other model parameters were fixed. This ensures 
the model structure is consistent with the assumed abundances. 



We determined T^s by fitting synthetic line profiles, calculated with SYNTH3 ([Kochukhov 



20071), to the observed profiles of two hydrogen lines: Hor (from the IDS spectrograph) and Hy 
(from the ESPaDOnS spectropolarimeter). We discarded the other hydrogen lines observed with 
ESPaDOnS because of the uncertainties in the continuum normalisation. In the temperature 
range expected for WASP- 12, hydrogen lines are extremely sensitive to temperature variations 



jhttp : //kurucz . harvard . edu 



and insensitive to logg variations, and ar e therefore 



Teff =6250+ 100 K, in good agreement with 



good t emperature indicators. We found 



HebbetaL 



(120091). The uncertainty estimate considered 



both the quality of the observations and the uncertainties in the normalisation. Figure [T] shows 
the comparison between the observed Har line profile and the synthetic profiles calculated with 
the adopted T^s of 625 0+ 100 K. The poor fit of the hydrogen line core is due to the adopted LTE 



approximation (see e.g. 



Mashonkina et al. 



2009) 



Another spectroscopic indicator for T^ff is given by the analysis of metallic lines, that we 
performed on the ESPaDOnS spectrum. In particular, T^ff is determined by eliminating the 
correlation between abundance and excitation potential O^fexcit) for the selected lines of a given 
ion/element. The top panel of Fig. |2] shows the correlation between abundance and Xexdt for all 
measured lines of Ca I, Ca II, Fe I, Fe II, and Ni I. Figure [2l produced using the final adopted 
fundamental parameters, shows no significant correlation between abundance and;t'excit for all 
ions, except Fe I, for which we registered a slightly positive correlation (0.02263+0.01 139), that 
would be eliminated by a higher Teff . What we found here for Fe I resambles what remarked by 



Ryabchikova et al. 



(120091) for HD 49933 and HD 32115 (both stars have a T^s slightly higher 
than WASP- 12): an effective temperature determination based only on the analysis of Fe I lines 
leads systematically to a T^ff that is substantially higher (by about 5%), compared to the one 
obtained with other ions and in particular with other temperature indicators, such as hydrogen 
lines U- For this reason, we decided to use the analysis of the metallic lines only as consistency 
check of Teff derived from the hydrogen lines. In addition to the ions shown in Fig.[2l we included 
in the consistency check also C I, Si I, Sc II, Ti I, Ti II, V I, Cr I, Cr II, Mn I, Co I, and Y II, 



^We do not know the precise origin of this phenomenon and also whether it is present in a 
large temperature range or just for F-type stars. We believe that more work should be done in this 
respect, in particular because the T^s determination based on the abundance-;^excit equilibrium is 
widely used to determine the eff"ective temperature of stars in a broad T^s range. 
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Fig. 1. — Comparison between the Hor line profile observed with the IDS spectrograph (black 
solid line) and synthetic profiles calculated with the final adopted Teff of 6250 K (red solid line), 
and uncertainty (dashed lines). 
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obtaining the requested equilibrium for all of them. The number of lines adopted to measure the 
abundance-;^excit correlation for each ion is the same as the one we used to derive the final ion 
abundance, and it is listed in Table [2J 

The surface gravity was derived from two independent methods based on (?) line profile 
fitting of gravity- sensitive metal lines wit h developed wings and ( //) ionisation balance for several 



elements. The first method is described in 



Fuhrmann et al. 



(|1997|) and uses the fact that the wings 



of the Mg I lines at /1 5167, 5172 and 5183 A are very sensitive to log^ variations. In practice we 
first derived the Mg abundance from other Mg I lines without developed wings, such as A 51 II 
and 5785 A, and then we fit the wings of the gravity indicator lines by tuning the logg value. To 
apply this method, very accurate log gf values and Van der Waals (log ywaai J damping constant s 



are req uired for all the lines. We adopted the set of line parameters used by 



Rvabchikova et al. 



(|2009|) and included the uncertainty in these pa rameters in t he unc ertainty in log^. We obtained a 



Hebb et al, 



(120091) . Our line profile fit of the Mg I 



logg value of 4.2+0.2, in good agreement with 
lines with developed wings is shown in Fig. [31 

The second method, the ionisation equilibrium, is often used to derive the surface gravity, 
but this method is extremely sensitive to the non-LTE eff'ects present for each ion/element, while 
Mg lines with developed wings ( less sensitive to non-LTE effects) are more suitable as logg 



indicators (Zhao & Gehren 



20001) . For this reason we decided to keep the Mg lines as our primary 



logg indicator, and checked the result with the ionisation equilibrium. Adopting the logg value 
obtained from the analysis of the Mg I lines with developed wings and taking into account the 
abundance uncertainties, we satisfied the ionisation equilibrium for every element with two 
analysed ions. 

Our main source for the atomic parameters of spectral lines is the VALD database with the 
default configuratio n. LTE abundan ce analysis was based on equivalent w idths, analysed with 



a modified version (Tsvmba 



1996 ) of the WIDTH9 code ( Kurucz 



1993a|) . For all analysed 
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Fig. 2. — Top panel: abundance obtained for each measured line of Ca I (open circles), Ca II (filled 
squares), Fe I (filled triangles), Fe II (stars), and Ni I (open downside triangles) as a function of 
the line;^fexcit- The black full line shows, as an example, the abundance-;^^excit correlation for Ni I: 
0.00637+0.01270, consistent with zero. Bottom panel: abundance obtained for each measured line 
of the ions given in the top panel, as a function of the measured line equivalent width. The elFective 
temperature from metallic lines is derived eliminating the correlation shown in the top panel, while 
f mic is determined eliminating the correlation shown in the bottom panel. 
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from 



Fuhrmann et al 



(l2009h: o scillator strengths from 



Aldenius et al. 



1997 ) 



(|2007h and damping constants 
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elements/ions we used almost all unblended spectral lines with accurate atomic parameters 
available in the wavelength range 4240-9900 A, except lines in spectral regions where the 
continuum normalisation was too uncertain. For blended lines, lines subjected to hyperfine 
splitting (hfs), lines situated in the wings of the hydrogen lines or for very shallow lines of 
specific ions we derived the line abundance by performing synthetic spectrum calc ulations with 



the SYN TH3 code. The hfs con stants for abundanc e calculations were taken from .Lawler et al 



(2001b) for Eu II lines and from 



Smith et al 



(119981) for the Li I line at A 6707 A. For the Barium 



abundance we use d the Ba II lines at /1 5853.7 A and 6496.9 A for which we do not expect any 



relevant hfs effect (IMashonkina & Zhao 



|2006() . A line-by-line abundance list with the equivalent 
width measurements, adopted oscillator strengths, and their sources is given in Table [T] (see the 
online material for the complete version of the table). 

The microturbulence velocity was determined by minimising the correlation between 
equivalent width and abundance for several ions, as shown in Fig. [2] for Ca I, Ca II, Fe I, Fe II, 
and Ni I. To employ this method we used all ions for which the measured spectral lines covered 
a large range in equivalent width. In particular we took into account simultaneously: Si I, Ti I, 
Ti II, Cr I, Cr II, Mn I, and Y II, in addition to the ones present in Fig. [2l The final adopted i/mic 
is 1 .2+0.3 km s~' . The given error bar is the range of values resulting from minimisation of the 
correlation for each ion we considered. 

The projected rotational velocity and macro turbulence (i/macro) were determined by fitting 
synthetic spectra of several carefully selected lines to t he ob served spectrum. Given the 



V sin i-Vrnacm dcgcncracy, we followed 



Valenti & Fished (|2005|) : we measured i/macro assuming 



f sin;=Okm s ^ and then vsini assu ming fniacro=4.75kms ^ (from the reff-i/macro relation 



published by 



Valenti & Fishei 



120051) . In the first case we obtained f macro =7. 0+0. 6 km s \ while 
in the second case we obtained f sin?=4.6±0.5kms '. In conclusion f sin? is in the range 



-1 



0-4.6+0.5 km s , while t'macro lays between 4.75 and 7.0+0.6 km s . Only with a careful analysis 
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Table 1: Linelist of the lines used for the abundance analysis. The first and second columns list 
respectively the ;t'excit (in eV) and the log gf value for each line. Columns five and six list the 
equivalent width in mA and the abundance for each line, while the last column gives the reference 
for the log gf value. Spectral lines for which the abundance was measured with synthetic spectra, 
instead of equivalent widths, present an "S" instead of the equivalent width value. Lines marked 
with "*" are subject to hyperfine structure, discussed in detail in the main text, while lines marked 
with "#" are multiplets (doublets or triplets) and in these cases we listed only the strongest line. 
The entire table can be viewed in the electronic version of the Journal. 
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of the Rossiter-McLaughlin (RM) effect will it be possible to precisely measure v sin /. 

The final WASP- 12 abundances, in \og{N/Ntot), are given in Table [21 a nd the atmospheric 



abundance pattern is shown in Fig.|4]in comparison to the solar abundances (|Asplund et al. 



20051) 



While the large ove rabundance of K would disappear if non-LTE effects are taken into account 



(Takedaet al. 



1996|) . the Sr overabundance is g enuine since non 



less then 0.05 dex for the Sr lines we analysed (|Mashonldna et al. 



TE ef fects are expected to be 



20071 and references therein). 



In general we expect small non-LTE effects for WASP- 12 due to the high stellar metallicity. 

The stellar metallicity (Z) is defined as follows: 

_ Z.>3m.lO'°g<^''/^-) 

where a is the atomic number of an element with atomic mass m^,. Our abundances imp ly a 



(1) 



Asplund et al. 



(2005) for all 



metallicity of Z=0. 02 1+0.002 dex, adopting the solar abundances by 
the elements that we did not analyse. 

The Z value adopted to characterise isochrones is calculated with the following 
approximation: 

^star — IVJ Z^O , 

assuming Zo=0.019dex. We recalculated the Z of WASP- 12 according to this approximation 
obtaining Z=0.036+0.002 dex. 



(2) 



4.1. Abundance uncertainties 

The abundance uncertainties for each ion, shown in Table [2l are the standar d deviations 
of the mean abundance obtained from the individual line abundances. Following 



Fossati et al. 



(120091) . it is possible to conclude that in case of ions with a sufficiently high number of lines, the 
internal scatter for each ion includes the uncertainties due to equivalent width measurement and 
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Table 2: LTE atmospheric abundances of WASP- 12 with error bar estimates based on the internal 



scatter from the number of analysed 
dex relative to the solar values from 
of the solar atmosphere from 



ines, n. The fourth co lumn gives the WASP- 12 abundances in 



Asplund et al. 



Asplund et al, 



(|2005h . The last column gives the abundances 



2005|) . The Lithium and Europium abundances take 
hyperfine structure in the lines into account. The Gd II abundance is an upper limit. The sym- 
bol # indicates the ions for which the abundance was derived from line profile fitting, instead of 
equivalent widths. 



Ion 


WASP- 1 2 




Sun 




logliV/Wioi) 


n 


[Wcl/Mal] 


loglN/Wioi) 


Lil# 


-9.47+0.05 


2 


+ 1.52 


-10.99 


CI 


-3.45+0.08 


14 


+0.20 


-3.65 


Nl# 


-3.95 


I 


+0.31 


-4.26 


01 


-3.10 


I 


+0.28 


-3.38 


Nal 


-5.63+0.04 


4 


+0.24 


-5.87 


Mgl 


-4.30+0.13 


4 


+0.21 


-4.51 


Mgll 


-4.29+0.05 


3 


+0.22 


-4.51 


All 


-5.63+0.08 


4 


+0.04 


-5.67 


Sil 


-4.47+0.18 


60 


+0.06 


-4.53 


Sill 


-4.33+0.01 


2 


+0.20 


-4.53 


Sl# 


-4.78+0.05 


8 


+0.12 


-4.90 


KI 


-6.12 


I 


+0.84 


-6.96 


Cal 


-5.48+0.07 


24 


+0.25 


-5.73 


Call 


-5.41+0.06 


7 


+0.32 


-5.73 


Sell 


-8.55+0.07 


8 


+0.44 


-8.99 


Til 


-6.95+0.09 


37 


+0.19 


-7.14 


Till 


-6.76+0.09 


24 


+0.38 


-7.14 


VI 


-7.99+0.07 


11 


+0.05 


-8.04 


VII 


-7.86+0.03 


2 


+0.18 


-8.04 


Cil 


-6.17+0.06 


32 


+0.23 


-6.40 


CiII 


-5.94+0.06 


12 


+0.46 


-6.40 


Mnl 


-6.41+0.16 


16 


+0.24 


-6.65 


Fd 


-4.31+0.12 


389 


+0.28 


-4.59 


Fell 


-4.22+0.06 


38 


+0.37 


-4.59 


Col 


-6.98+0.07 


12 


+0.14 


-7.12 


Nil 


-5.64+0.10 


105 


+0.17 


-5.81 


Cul 


-7.81+0.09 


4 


+0.02 


-7.83 


Znl 


-7.32 


I 


+0.12 


-7.44 


Srl# 


-8.30 


I 


+0.82 


-9.12 


SrII# 


-8.35+0.05 


3 


+0.77 


-9.12 


YII 


-9.55+0.09 


10 


+0.28 


-9.83 


Zrll 


-9.08 


I 


+0.37 


-9.45 


Ban 


-9.37+0.06 


2 


+0.50 


-9.87 


Lall 


-10.40 


I 


+0.51 


-10.91 


Cell 


-10.29+0.07 


3 


+0.17 


-10.46 


NdU 


-10.45+0.03 


4 


+0.14 


-10.59 


Smll 


-10.64 


I 


+0.39 


-11.03 


EuII# 


-11.30 


I 


+0.22 


-11.52 


GdU# 


< -10.92 


I 


+0.00 


-10.92 


J-eff 




6250 K 




5777 K 


logu? 




4.20 




4.44 



18 



X 

3 



> 

-(— > 

u 

O 

a 



1.6 
1.4 




























































( 


t 
























































- 


- 




























































- 


1.2 


• neutrals 

A singly ionized 




- 














































- 


1 

0.8 
0.6 
0.4 
0.2 



















- 
















« 


» 






















\ 


► 
















- 


- 








































L 


i 
















- 


- 






















L 


^ 


- 


Z 




L 
\ 














- L 


L 


-. ^ 


\ 


[ LA 


_ 


- 




( 




» 


[] 


k.... 
( 


■ z 






i 
\ 


r 


_ L 


i 


( 


» 


- 


- 






I 




' 


- 






( 

» - 


» 






< 


1 L 


li - 

1 




- 


« 


( 


» 
( 


_ « 


» 










L\ -- 




/. 


- 














- 
























- 


















V 



element 



Fig. 4. — Ion abundance relative to the Sun (lAsplund et al. 



2005h of the WASP- 12 atmosphere. Full 



points show the abundance of the neutral elements, while the open triangles show the abundance 
of the singly ionised elements. Each abundance value is shown with two uncertainty values: the 
standard deviation from the mean (column 3 of Table [3]) and the total uncertainty (column eight of 
Table O. We believe the real uncertainty lies between these two values. 
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continuum normalisation. In addition, from plotting the abundance scatter as a function of the 
number of lines, we can also infer an internal uncertainty of 0.08 dex, that can be applied as mean 
scatter when only one line of a certain ion is measured. 

Table [3] shows the variation in abundance for each analysed ion, caused by the change of one 
fundamental parameter by +lcr, keeping fixed the other parameters. 

Table [3] shows that the main source of uncertainty varies according to the element/ion (e.g. 
for the Fe-peak elements, neutrals are more sensitive to temperature variations, while ions are 
more sensitive to log^ variations) and in some cases to the selected lines (e.g. the two Ba II lines 
selected to measure the Ba abundance are rather strong, therefore the Ba abundance is strongly 
dependent on Vmic variations). 

Assuming the different uncertainties in the abundance determination are independent (though 
actually the systematic uncertainties will be correlated), we derived a pessimistic final error bar 
using standard error propagation theory, given in columns seven and eight of Table (3] Using the 
propagation theory we considered the situation where the determination of each fundamental 
parameter is an independent process. The mean value of the LTE uncertainties given in column 
eight of Table [3] is 0.1 1 dex. Due to the fact that for the parameter determination of both T^s, 
logg and fmic we took into account all possible systematics (except non-LTE), we believe that the 
abundance uncertainties given in the last column of Table [3] can be considered as upper limits and 
that the real error bars lie between the values given in column three and eight of Tabled 

5. High precision magnetic field search 

One of the main goals of our analysis is to search for a global stellar magnetic field in 
WASP- 12. The ESPaDOnS spectropolarimeter yields high resolution and high SNR spectra 
of both Stokes / and V allowing this search. To detect the presence of a global magnetic field 
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Table 3: Uncertainty sources for the abundances of WASP-12. The third column shows the stan- 
dard deviation (Xabn (scatt.) of the mean abundance obtained from different spectral lines (internal 
scattering); a blank means that only a single line was used, and we estimated the internal scattering 
to be 0.08 dex. Note that these values are identical to those given in Table [2l Columns 4, 5, and 6 
give the variation in abundance estimated by increasing Teff by 100 K, log g by 0.2 dex, and i/mk by 
0.3kms \ respectively. Column 7 gives the the mean error bar calculated adding the systematic 
uncertainties given in columns 4, 5 and 6 in quadrature i.e., cr^^^^ (syst.) = crl^^ (Tsff) + crl^^ (log^) 
+ crl^^ (f mic)- Column 8 gives the total mean error bar: o-^^^ (tot.) = o-^^^ (syst.) + (rl^^^ (scatt.). 





Ion 


abundance 


iTabn (scatt.) 


o-abn (Ten) 


D-i,bn(logS) 


0"abn ("mit) 


o-abn (sysl.) 


o-abi, (lot.) 




logfW/Mol) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


Lil 


-9.47 


0.05 


0.06 


-0.01 


0.00 


0.06 


0.08 


a 


-3.45 


0.08 


-0.05 


0.05 


0.00 


0.07 


0.11 


NI 


-3.95 




-0.06 


0.05 


0.00 


0.08 


0.11 


01 


-3.10 




-0.08 


0.05 


0.00 


0.09 


0.12 


NaT 


-5.63 


0.04 


0.04 


-0.04 


-0.02 


0.06 


0.07 


Mgl 


-4.30 


0.13 


0.04 


-0.02 


-0.02 


0.05 


0.14 


Mgll 


-4.29 


0.05 


-0.08 


0.04 


-0.02 


0.09 


0.10 


All 


-5.63 


0.08 


0.03 


-0.01 


-0.01 


0.03 


0.09 


Sil 


-4.47 


0.18 


0.02 


-0.01 


-0.01 


0.02 


0.18 


Sin 


-4.33 


0.01 


-0.07 


0.06 


-0.02 


0.09 


0.09 


SI 


-4.78 


0.05 


0.06 


-0.06 


0.00 


0.08 


0.10 


KI 


-6.12 




0.08 


-0.08 


-0.04 


0.12 


0.14 


Cal 


-5.48 


0.07 


0.07 


-0.05 


-0.06 


0.10 


0.13 


Can 


-541 


0.06 


-0.06 


0.04 


-0.02 


0.07 


0.10 


Sen 


-8.55 


0.07 


0.01 


0.07 


-0.04 


0.08 


0.11 


TU 


-6.95 


0.09 


0.08 


0.00 


-0.02 


0.08 


0.12 


■nn 


-6.76 


0.09 


0.00 


0.07 


-0.06 


0.09 


0.13 


VI 


-7.99 


0.07 


0.09 


0.00 


-0.01 


0.09 


0.11 


vn 


-7.86 


0.03 


0.00 


0.08 


-0.01 


0.08 


0.09 


QI 


-6.17 


0.06 


0.06 


-0.01 


-0.02 


0.06 


0.09 


Crn 


-5.94 


0.06 


-0.02 


0.08 


-0.03 


0.09 


0.11 


Mnl 


-641 


0.16 


0.07 


0.00 


-0.03 


0.08 


0.18 


Fel 


-4.31 


0.12 


0.06 


-0.02 


-0.04 


0.07 


0.14 


Fen 


-4.22 


0.06 


-0.02 


0.07 


-0.04 


0.08 


0.10 


Col 


-6.98 


0.07 


0.07 


0.00 


-0.01 


0.07 


0.10 


Nil 


-5.64 


0.10 


0.06 


-0.01 


-0.03 


0.07 


0.12 


Cul 


-7.81 


0.09 


0.06 


0.00 


-0.02 


0.06 


0.11 


Znl 


-7.32 




0.02 


0.02 


-0.02 


0.03 


0.09 


Sri 


-8.30 




0.02 


-0.01 


-0.03 


0.04 


0.09 


Srn 


-8.35 


0.05 


0.00 


0.02 


-0.02 


0.03 


0.06 


Yn 


-9.55 


0.09 


0.01 


0.08 


-0.03 


0.09 


0.12 


an 


-9.08 




0.01 


0.08 


-0.01 


0.08 


0.11 


Ban 


-9.37 


0.06 


0.04 


0.02 


-0.11 


0.12 


0.13 


Lan 


-10.40 




0.03 


0.08 


-0.01 


0.09 


0.12 


Cen 


-10.29 


0.07 


0.03 


0.08 


-0.01 


0.09 


0.11 


Ndn 


-10.45 


0.03 


0.03 


0.08 


-0.01 


0.09 


0.09 


Smn 


-10.64 




0.02 


0.08 


-0.02 


0.08 


0.12 


EuII 


-11.30 




0.03 


0.08 


-0.01 


0.09 


0.12 


Gdn 


< -10.92 




0.02 


0.08 


0.00 


0.08 


0.11 
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and measure its strength we used the Least-Squares Deconvolution technique (hereafter LSD), 
adopting a code written by one of us (O. Kochukhov). 

LSD is a cross-correlation technique developed for the detection and measure ment of weak 



polaris atio n signatures in stel lar spectral lines. The method is described in detail by 



( 1997 ) and 



Wade et al. 



M- 



Donati et al 



(|2000l) . We decided to use the LSD approach to detect a magnetic field in 
WASP- 12 since this method is the most precise currently available, especially for stars with rich 
line spectra and low projected rotational velocity (vsini), such as WASP- 12. 

We applied the LSD technique to the Stokes V spectra from each of the two nights using 
about 6 450 atomic spectral lines with the only cut-off" criterion based on the calculated line depth 
(>10%), and in both cases no magnetic field was found. From the spectrum of January 3rd we 
obtained (B,) = 2.3 + 5.3 G, with a SNR of the Stokes V LSD profile of 2400, while from the 
second night spectrum we obtained (B^) = -10.1 ± 4.2 G, with a SNR of the Stokes V LSD profile 
of 3380. We obtained similar values also from both null profiles. Figure |5] shows the LSD profiles. 



In late-type stars, the stellar magnetic field is directly connected to the chromospheric 
activity, that can be monitored with the Ca II H and K lines. In Fig. [6] we compare the profiles 
of the Ca II H and K lines, observed with ESPaDOnS, with the mean line profiles of tBoo 



(Shkolniketal. 



2005|) obtained averaging several CFHT spectra acquired with the GECKO 



spectrographo. This comparison is particularly valuable because both stars are planet hosting and 



have similar fundamental parameters and age (JGonzalez et al. 



2010a[) . where the difference is 



mainly in the v sin i values (r Boo has a v sin i of about 13.5 km s ^). Figure [6] does n ot show the 



presence of any anomaly in WASP- 12 Ca II H and K line cores. 



Knutson et al. 



20101) determined 



the \og{R'r,j^) chromospheric stellar acitvity parameter in a set of planet hosting stars, reporting for 



'\http://w ww. cflit. hawaii. eduJInstruments/Spectroscopy/Gecko^ 
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Fig. 5. — LSD profiles from tlie spectra obtained on January 3rd (top panel), corresponding to an 
orbital phase of 0.80, and January 5th 2010 (bottom panel), corresponding to an orbital phase of 
0.81. In each panel, the top black line shows Stokes /, normalised to the continuum (the corre- 
sponding units are given in the right y-axi s). The middle red line corresponds to the V/I profile 



and the bottom blue line to the null profile (IDonati et al. 



19971); their corresponding units are given 



in the left y-axis. The Stokes / LSD profile is centered on the stellar radial velocity and the null 
profile is shifted downward by an arbitrary offset. 
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WASP-12 the remarkably low value of -5.500, the lowest in their sample. 



The bottom panel of Figure 1 in 



Fossati et al. 



(120101) shows the core of the Mg II UV 



resonance lines (further stellar activity indicators) of WASP-12 and here the lac k of any line co re 



emission is clear, in agreement with the low level of stellar activity reported by 



Knutson et al 



( 2010 ) 



It is well known that the chromospheric activity is strongly correlated with the stellar 



rotational ve 



Husnoo et al. 



ocity, t hat for WASP-12 is unknown. Given the very small RM effect shown by 
(|2010l) . it is probable that the low stellar activity of WASP-12 might be connected 



to a low rotational velocity. It is also possible that the sta r is either passing through a period of 



low activity, or WASP-12 is much older than reported by 



Hebb et al. 



(l2009h (more than 1 Gyr 



older), or the typical signs of stellar activity, such as emission in the cores of the Mg II resonanc e 



lines, are absorbed by the material lost by the planet and falling onto the star (see 



Fossati et al 



m 



Li et al 



20101 : 



20101) . A more thorough analysis and discussion of the stellar activity and of the 



reasons behind its low level will be given in a following work, now in preparation. 

6. Discussion 

6.1. Atmospheric parameters and convection 

Because WASP-12 has a relatively low effective temperature, its hydrogen lines are ideal 
indicators of atmospheric T-P structure. Generally, fitting hydrogen line profiles, rather than any 
other atomic line, provides an accurate estimate of T^fi. On the other hand, at low temperatures, 
convection becomes an important energy transfer mechanism, influencing the photospheric 
temperature stratification and, thus, hydrogen line formation. The derived value of T^ff then 
depends on the convection treatment used in the model atmosphere calculations. 



There are basically two formalisms of convection treatment available for model atmosphere 
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Fig. 6. — Comparison of the observed Ca K and H line (left and right panels ^ respectively) profile s 



Shkolnik et al. 



2005 ). 



between WASP- 12 (thin line) and the planet hosting star tBoo (thick line, 

The two upper panels show the WASP- 12 spectrum obtained on January 3rd, 2010, while the two 

lower panels show the spectrum obtained on January 5th, 2010. 
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computations: the well known mixing-length theory (MLT) (lBohm-Vitenselll958|) . which relies on 
some free parameter (or), describing the characteristic length traveled by co nvective cells before 



they d isappear, and further improvement of stellar convection developed by 



Canuto & Mazzitelli 



(|l99ll CM hereafter). The main advantage of the new model is that it does not require any 



adjustable parameters like a, which is now computed based on the geometrical depth scale 
inside the convective zone. It also accounts for the full spectrum of turbulent eddies, and is 
therefore superior to the sin gle-eddy assumption made by MLT. More details can be found in 



Canuto & Mazzitellil (119921) . 



Physically the CM convection model is somewhat superior to MLT, therefore we preferred 
CM convection over MLT in our analysis. In the case of WASP- 12, we find no critical differences 
in Ha profiles using both CM and MLT convection theories with the commonly accepted a = 1 .25. 
This is because the wings of Ha are formed mainly in the region right above the photosphere 
where convection plays a less important role in the energy balance. However, bluer Balmer lines 
are formed deeper and thus can show the changes in the temperature stratification introduced by 
convection. For instance, this is the case for the Hy line for which MLT predicts weaker line 
wings (i.e. stronger convective energy transport): one needs an approximately 100 K hotter model 
to fit the Hy line wings with MLT convection. However, it is then clear that only CM allows a 
consistent fit simultaneously for Ha and Hy lines with the same T^g. 



6.2. The distance and age of WASP-12 



The basic parameters of WASP-12 were derived by .Hebb et al. 



mm. Teff, log^. 



metallicity, and v sin i were derived from the observed stellar spectrum using spectral synthesis; 
a simultaneous Markov chain Monte Carlo (MCMC) fit to radial velocity and transit light 
curve measurements produced values for Rp/R^, R*, M* and the orbital semi-major axis a. The 
observed spectral type of the star and its radius imply the luminosity and the spectral energy 
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distribu tion. Using this inform ation, the 2MASS IR fluxes can be used to derive the distance 



(IRFM; 



Blackwell et al. 



19791). This method yields a distance to WASP-12 of 265+20pc 



(|http : //www ■ superwasp . org/wasp_planets . htm|), assuming a typical main-sequence ste llar 



Hebb et al. 



( 2009 ) 



radius; and 385+30 pc, assuming the stellar radius from the MCMC fits of[ 
(B. Smalley, priv.comm.). We note that the IRFM implicitly assumes solar metallicity and zero 
reddening, both inapplicable for WASP-12. Since the orbital period and the radial velocity 
amplitude fix the value of a, a nd the trans i t light curve fixes the ratio a:i?,, the stellar distance of 



385+30 pc is to be preferred. 



Hebb et al 



(|2009|) also estimated the age of WASP-12 by fitting 



isochrones to the position of WASP-12 in a modified HR diagram (temperature vs. inverse cube 
root of the stellar density), obtaining an age of 2.0;^^^ Gyr old; given the several uncertainties, the) 
increased the error bars to 1 Gyr, concluding that WASP-12 is 2+1 Gyr old. 



One of the most secure empirical properties of the star is its colour, or equivalently effective 
temperature. We can assess the distance and age of WASP-12 by comparing the eff"ective 
temperature with is ochrones and evolutionary tracks. Figure |7] shows four isochrones from 



Marigo et al. 



(|2008l) with a metallicity Z of 0.030, the maximum available value of Z, adopted 
following Eq. |2l WASP- 12's eff"ective temperature places it on the vertical thick blue line of Fig. Ul 
where the other two full vertical blue lines are defined by the uncertainty on T^s- On this central 
line we placed three points, which correspond to distances of 265 pc (lower point), 295 pc (middle 
point) and 465 pc (upper point). The lower point lies well below the zero-age main sequence 
(ZAMS). From this we can rule out a distance as close as 265 pc; this concurs with the MCMC 



Hebb et al. 



(|2009|) in implying the star is bigger than the typical main sequence stellar 



fitting of 
radius. 

We can see that a range of possible ages and distances are compatible with the isochrones. 
In the region of Fig. |7]between the two solid triangles, several isochrones are consistent with the 
empirical eff"ective temperature. The 2.65 Gyr isochrone is just consistent with the lower limit 
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Fig. 7. — Position of WASP-12 on the HR diagram assuming three different stellar distances: 
265 pc (circle), 295 pc (triangle), and 465 pc (inverted triangle). The maximum and minimum dis 



tances were calculated adopting Tgff =6250 K and interstellar reddenin gs from (lAmores & Lepine 



20051). The dotted, thin full and dashed lines show isochrones from 



Marigo et al. 



(l2008h corre- 



sponding to ag es of 1 Gyr, 2 Gyr and 3 Gyr, respectively, encompassing the possible age range of 



WASP-12 from 



Hebb et al. 



(|2009|) . The thick full line is the 2.65 Gyr isochrone we a rgue this is the 



Girardi et al 



maxim um possible age for WASP-12. The red lines show evolutionary tracks from 
hood) for 1.5M/Mq, 1AM/Mq, 1.3M/Moand 1.2 M/Mq, from top to bottom. Both isochrones 
and evolutionary tracks assume a metallicity Z of 0.03. The blue vertical lines show the WASP- 
12's temperature range; these lines change from full to dashed below the ZAMS, indicated by the 
green line. 
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on WASP-12's effective temperature: the loop to the right at the main sequence tum-ofF point of 
this isochrone intersects the lower limit on the temperature. No star of this age is consistent with 
WASP-12's colour except for those which are turning on to the horizontal branch, higher up in 
the diagram at logL/Lo«l. For stars at this stage of evolution log^«3.8, whereas the spectrum of 
WASP-12 implies logg~4.2 (c.f. Sect|4]). WASP-12 thus cannot be turning on to the horizontal 
branch. Consequently we can constrain the position of WASP-12 in Fig.|7]to be around or younger 
than the main sequence turn-off. 

Applying this reasoning, the oldest possible age for WASP-12 is 2.65 Gyr. This arises from 
the intersection of lower limit on WASP-12's effective temperature and the full thick isochrone 
in Fig. |71 For all isochrones younger than this, there is also an intersection at or before the 
main-sequence turn-off. For the isochrones older than this no allowed intersection occurs: only 
evolved stars have compatible temperatures, but these are ruled out by their surface gravity. 



Girardi et al 



(|2000l) . for stars of mass 



The full red lines in Fig. |7] are evolutionary tracks, from|j 
1.2M/Mo, 1.3M/Mo, 1.4 M/Mq and 1.5 M/Mq respectively from bottom to top. WASP-12 is 
clearly hotter than a 1 .2 M/Mq star for any age, therefore we conclude its mass exceeds 1 .2 M/Mq. 
Interpolating between the evolutionary tracks, we estimate a limit on the mass of WASP-12 of 



(120091) . but is a tighter 



around 1.23 M/Mq. This is consistent with MCMC fitting of .Hebb et al. 
constraint. 

The evolutionary track for 1.3 M/Mq almost exactly coincides with the intersection of the 
2.65 Gyr isochrone with the lower limit on the effective temperature. The oldest possible age 
therefore corresponds to a stellar mass of about 1.3 M/Mq. For higher masses, the empirical 
effective temperature intersects the evolutionary track while the star is still on the main sequence, 
and is therefore consistent with WASP-12's logg. The 1.4 M/Mq evolutionary track is consistent 
with the empirica l effective temperature, and is also consistent with the MCMC fitting of 



Hebb et al 



(l2009h . 
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The upper limit on the mass from lHebb et alj (|2009h is just below 1.5 M/Mq. We can use 
this to infer an upper limit on the luminosity, and hence the distance. We find that the maximum 
distance is 465 pc. On the other hand, the minimum distance is obtained at the intersection of 
the vertical line that defines WASP-12's effective temperature and the ZAMS; this distance is 
295 PC. These distances are computed taking into account insterstellar reddenings as given by 



(|Amdres & Lepind 



2005 ). 



One effective way to measure the age of a late-type star is the comparison of the lithium 
abundance with that of open cluster member stars for which the age is precisely known. For a 
more thorough analysis of the Li I line at A 6707 A we downloaded from the SOPHIE archivqj all 
the 21 mid-r esolution spectrcQ (R~40000) of WASP- 12 obtained for the radial velocity analysis 



published by 



Hebb et al. 



( 200' 



i. 



Figure [8] shows a comparison between the 21 SOPHIE spectra and synthetic spectra 
calculated with the adopted stellar parameters and abundances around the region of the Li I line at 
A 6707 A; the bottom-right panel shows the same comparison, but with the ESPaDOnS spectrum. 
This plot shows that the Li abundance derived from the ESPaDOnS spectrum fits the SOPHIE 
data as well and that there is no line profile variation of the Li line with the orbital phase. We 
measured also the equivalent width of the Li line in each of these spectra and we did not find any 
significant time variation. For this reason we believe that t he lac k of detection of this Li I line in 



Hebb et al 



(120091) ■ could be due to the low SNR or 



the SARG spectrum of WASP- 12, reported by 
to a wrong line identification. 

From a comparison between the lithium abundance of WASP- 12 and the results published 



^httpV/atlas.obs-hp.fr/sophie/ 

^SOPHIE is a cross-dispersed echelle spectrograph mounted at the 1.93-m telescope at the 
Observatoire de Haute-Provence (OHP). 
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Fig. 8. — Comparison between the observed phase dependent SOPHIE spectra (black solid line) 

o 

of WASP- 12 around the Li I line at /1 6707 A and synthetic spectra (dashed red line) calculated 
with the Li abundance and the parameters obtained from the analysis of the ESPaDOnS data. 
The bottom-right panel shows the same comparison, but with the analysed ESPaDOnS spectrum 
instead. The synthetic spectra take into account the difference in resolution between SOPHIE and 
ESPaDOnS. 
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by lSestito & RandichI (|2005|) we can only conclude that WASP- 12 is older than 500 



is in agreement with the 2.0+1.0 Gyr given by evolutionary tracks (|Hebb et al 



Myr, which 



20091) . The large 



uncertainty is caused by the extremely slow lithium depletion for stars with effective temperatures 
similar to that of WASP- 12, due to the shallow surface convection zone. 

In conclusion o ur analysis on the age of WASP-12 leads to a stellar age between 1.0 Gyr 



(lower limit given by 



Hebb et al. 



2009|) and 2.65 Gyr (from the analysis of the HR diagram), as 



shown also in Table |4] that lists the minimum and maximum values we obtained for the stellar age, 
mass and distance. 

Table 4: Minimum and maximum values of WASP-12's age (in Gyr), M/Mq and distance (in pc) 
as derived by the analysis of the HR diagram. The values marked with an "*" are taken from 



Hebb et al. 



(l2009h . 



Parameter 


Minimum 


Maximum 




value 


value 


age (Gyr) 


1.0* 


2.65 


M/Me 


1.23 


1.49* 


distance (pc) 


295 


465 



6.3. Is WASP-12 a chemically peculiar star? 



One of the main purpose of this work is to search for chemical peculiarities that could be 
connected with pollution of the stellar atmosphere by material lost by the planet. T his can be do ne 



in dif ferent ways, as also shown in the extensive salient literature present (see e.g. 



Neves et al. 



20091 and references therein). WASP-12 is a promis ing target for signs of atmospheric pollution: 



WASP-12 b is most likely currently losing material (|Fossati et al. 



2010|) : this material is believed 
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to be forming a circumstellar disk that is accreting onto the star, polluting the stellar photosphere 



(ILi et al 



2010|) . Classical ways of looking for atmospheri c pollution are 



peculiarities of single elements, such as Li and Be (e.g. 



Israelian et al. 



?V sea rching: (?) chemical 



the element abundance against the condens ation temperature (T^; see e.g. 



20041): (//) a trend in 



Sadakane et al. 



Ecuvillon et al. 



2006 



Melendez et al. 



20021 : 



20091) : (///) chemical peculiarities of the abundance pattern 



in comparison to reference stars (both planet hosting and non-planet hosting). 

When searching for small effects on the stellar atmospheric abundances it is important to 
determine whether a certain star belongs to the thin or thick Galactic disk population. To do so 
we performed both a kinematic and chemical analysis of WASP- 12. We calculated the Galactic 
velocity vectors (U, V^ W) co rrected to the local standard of rest (LSR) using the formalism of 
Johnson & SoderblomI (Il987h. in stead defining U as positive towards the Galactic anti-centre. 



As done by 



Sozzetti et al 



Soubiran & Girard 



(|2006r) . WASP- 12 was then placed on the Toomre diagram of the 



(|2005|) stellar sample, indicating that WASP- 12 has a peculiar velocity less 



than 85 km s , strongly indicative of thin disk membership. From the WASP- 12 abundances, we 



then compared our values of [a/Fe] (-0.15 dex) and 
as 0.25([Mg/Fe]+[Si/Fe]+[Ca/Fe]+[Ti/Fe]), with the 



Fe/H] (-1-0.28 dex)^ where [a/Fe] is defined 



Soubiran & Girard 



(|2005h sample, obtaining 



that WASP- 12 is again consistent with the properties of the thin disk population. 

In our analysis we are not able to look for chemical peculiarities of both Li and Be. As 
shown in Fig. [81 our spectra do not have the SNR and the resolution necessary to perform a precise 
analysis of the Li^/Li^ ratio, but we c an compare the lithiu m ab undance obtained for W ASP- 12 



with that of other planet hosting stars. 



Israelian et al. 



(2004) and 



Gonzalez et al. 



(l2010al) published 



lithium abundances of stars hosting a giant planet or a brown dwarf and compared them to a set of 
non-planet hosting stars. The lithium abundance we obtained for WASP- 12 matches that of both 
set of stars, showing clearly that Li is not peculiar in WASP- 12. Our spectra do not include the 
region around /1 3 1 30 A covering the two Be II lines usually adopted to measure the Be abundance. 
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6.3.1. Volatile vs. refractory elements 

An effective way to check whether a stellar photosphere is polluted by accretion of metal-rich 
material is to examine the correlation between the relative abundance of various elements and 
their T^. In accreting metal-rich material the refractory elements tend to form dust grains that are 



blown away by the stellar wind. Thus the s tar accretes more volati 



as happens, for example, in A Bootis stars (IMartinez-Galarza et al 



e than refractory elements. 



2009 LOur WASP- 12 ion 



abundance are shown against condensation temperature (taken from .Lodders 



20031) in Fig. m 



(e.g. 



The correlation obtained between the ion abundances relative to the Sun (lAsplund et al. 
20051) and T^ is not statist i cally significant: 6.91+5.64x10'^. In the past years several authors 



Sadakane et al. 



20021 : 



Ecuvillon et al. 



2006 



Sozzetti et al. 



20061) looked for the same kind 



of correlation in several planet hosting stars, without having found any. The trend we obtain for 
WASP- 12 is also indicative of a null result. 



Melendez et al.l (|2009|) . applying a differential analysis on the Sun and solar twins, discovered 
that the Sun shows a highly significant correlation of the abundances (relative to the mean 
abundance of the solar twins) in respect to the con densation temperature. In a ddition this 
correlat ion shows a significant break at T^ ~1200 K. Melendez et al.l (|2009[ ) and lGonzalez et i 
(l2010b|) came to the conclusion that this correlation is likely connected to the presence of planets 
around the Sun. The detection of this kind of correlation requires very high precision abundances 
(error bars < 0.05 dex and free from systematic differences), attainable only with a differential 
analysis. The application of this method to WASP- 12 would require observations of stars that can 
be considered twins of WASP- 12 (similar Teff, metallicity, age, and population). The choice of the 
WASP- 12 twins will be of crucial importance because effects such as diffusion, more prevalent for 
F-type stars than G-type stars, are strongly age and Tes dependent, and consequently could easily 
hide or mimic pollution signatures. 



Our analysis does not lead to a firm conclusion about the atmospheric pollution of WASP- 12. 
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Fig. 9. — Ion abundance relative to the sun (,Asplund et al. 



2005 ) as a function of the condensation 



temperature. The full line shows the linear fit to the data. The abundance uncertainties are given 
as in Fig. m 
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6.3.2. Comparison with previous analysis 

To check whether WASP- 12 has a peculiar abundance pattern we compared the abundances 
with those obtained by other authors on a large number of stars with a similar Teff. In particular 
we took into account the results of the abundance analysis on stars that are not known planet 
hosts, in the temperature range 6000-6500 K, and with a [Fe/H]>O.Odex. We decided to use 
such a small temperature range to decrease the effect of possible systematics, such as non-LTE 
effects. In addition, since we are interested in looking for a peculiar pattern, we decided to use 
only the stars with an over-solar iron abundance and to subdivide the sample of comparison stars 
according to their iron abundance: 0.0<[Fe/H]<0.1, 0.1<[Fe/H]<0.2, and [Fe/H]>0.2. To have 
a better statistical view of each subsample we decided to plot in our comparisons not only the 
mean abundance and relative standard deviation for the comparison stars, but also their abundance 
range. To accomplish all these re quirements, we needed t o have a large set of non-planet hosting 



stars, such as the one provided by 



Valenti & Fished (|2005h 



Valenti & Fisheij (|2005|) performed a parameter determination and LTE abundance analysis 
of more than 1000 late-type main sequence stars. For each star they derived the abundances of 
Na, Si, Ti, Fe, and Ni. This sample provided us with comparison abundances for 59 stars with 



0.0<[Fe/H]<0.1 dex, 36 stars with 0.1<[Fe/H]<0.2dex, and 22 stars with [Fe/H 
compa rison between our WASP- 12 abundances and the ones from the sample of 



>0.2dex. The 



Valenti & Fisher 



jlOOSh is shown in Fig. [JOl 



Figure \T0\ shows that the abundances of Ti and Fe are clearly comparable to the ones of the 
higher metallicity subsample, as expected given the fact that [Fe/H] of WASP- 12 is about 0.3. 
The Na abundance of WASP- 12 falls within the abundance range of the high-metallicity stars of 
the comparison sample, due to the large range covered by the Na abundance, while Ni shows a 
mean abundance that is just outside the abundance range. The Si abundance looks definitely more 
similar to the one of the set of stars with 0.1<[Fe/H]<0.2 dex, but the large uncertainty on this 
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Fig. 10. — Comparison between the mea n abundance, relat i ve to t he sun, of WASP-12 (full circles) 



and of the selected sample of stars from 



Valenti & Fished (|2005|) . that was subdivided according 



to the iron abundance: 0.0<[Fe/H]<0.1 dex (open upside triangles), 0.1<[Fe/H]<0.2dex (open 
squares), and [Fe/H]>0.2dex (open downside triangles). The shaded areas show the abundance 
range, while the error bars of the comparison sample give the standard deviation from the mean 
abundance. The uncertainties on the abundances of WASP-12 are as in Fig.lH 



37 



abundance does not allow any firm conclusion. In summary the abundances of Ti and Fe follow 
the pattern of the high-metallicity stars, while Na, Si, and Ni seem to follow more the abundance 
pattern of the set of comparison stars with [Fe/H] between 0.1 and 0.2 dex. 

The information gathered from this comparison lead to the conclusion that there is the 
possibility that the WASP- 12 abundances of Na, Si, and Ni do not really match the abundance 
pattern that the WASP- 12 iron abundance would suggest. But it has to be taken into account 
that the comparison is not free from systematic differences in the fundamental parameters and 
abundance determinations, and in the set of adopted atomic line parameters. 

For a further search of possible abundance peculiarities in the atm osphere of WASP 



we decided to compare our results w'ith those of several other authors (ITakeda et al 



Sadakane et al. 



2006 



2002; 



Bond et al 



Ecuvillon et a 



20081 : 



Neves et al 



20M; 

T 



Sozzetti et al. 



20061 : 



Gilli et al. 



2006 



2002 



Santos et al 



2009() have obtained in planet hosting and non-planet hosting 



stars. In this way we should also be able to partially remove possible systematic differences in 
case a peculiar pattern becomes evident in the comparisons with several different authors. 

The WASP- 12 abundances we derived for O, Ca, Sc, Mn, Zn, Y, Zr, Nd, and Eu match very 
well the ones previously obtained by other authors in planet hosting and non-planet hosting stars, 
while abundances of Al, V, and Cu appear to be clearly below the values previously obtained by 
more than 0.2 dex. For the other compared elements (C, N, Na, Mg, Si, S, Ti, Cr, Co, and Ni) we 
obtained just a satisfactory agreement, since we register a tendency of the WASP- 12 abundances 
of these elements to lay always in the lower margin of the comparison samples, in particular when 
the abundances are put in relation to the iron abundance. This result follows what previously 
obtained, strengthening the possibility of an increased Fe abundance in comparison to the one of 
some other elements. 



by 



In addition we compared our Si, Ti, and Ni WASP- 12 abundances with the ones obtained 
in a set of planet hosting and non-planet hosting stars. We obtained a 



Robinson et al 



(I2006D 
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good agreement for Ti, but both Si and Ni app ear to be depleted in WASP- 12, compared to the 



abundances obtained by 



Robinson et al 



(|2006h in a set of planet hosting stars, therefore we cannot 



confirm their conclusion of a systematic enrichment of Si and Ni in planet hosting compared to 
non-planet hosting stars. 



6.4. Spectral energy distribution: searching for a circumstellar disk 



The circumstellar disk, predicted by lLi et al.l(|2010h and tentatively observed by iFossati et al 



(|2010l) . could be detectable in the infrared. For this reason we decided to compare the 



calibrated Near-UV COS fluxes of WASP- 12 and the available photometry (Johnson and 2MASS 
photometry) with synthetic fluxes obtained with LLmodels adopting the fundamental parameters 
and the abundances derived from the ESPaDOnS data, looking for an infrared excess in the 
2MASS photometry. For this comparison, sh own in Fig. [TTJ we took into account a reddening 



E(B-V)=0.126mag (|Amdres & Lepine 
obtained in Sect. 



20051), that was calculated assuming the stellar distance 



Figure \TT\ shows a very good agreement between the model fluxes and the observations, 
demonstrating also the quality of both the adopted fundamental parameters and interstellar 
reddening. 

Figure [TT] does not show any clear infrared excess in the 2MASS photometry that could be 
interpreted as: (z) there is no interstellar disk around the star or (ii) the interstellar disk is not 
visible becau se either t he em ission is not strong enough or the emission peaks at much longer 



wavelengths. 



Li et al. 



(|2010|) suggested that the disk emission should peak at ~2.3//m, inside 
the wavelength range covered by the K band 2MASS photometry. Therefore, to look for an 
interstella r disk around W ASP- 12 it would be necessary to observe the system with a very high 



precision (|Li et al. 



2010h and at longer wavelengths, both accessible with the Herschel satellite 
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Fig. 11. — Comparison between LLmodels theoretical fluxes calculated with the fundamen- 
tal parameters and abundances derived for WASP- 12, taking into account a reddening of E(B- 
V)=0.126mag (full line) and without taking into account reddening (dashed line), with COS cali- 
brated fluxes (thick gray line), Johnson BV photometry (full points) and 2MASS photometry (full 
triangles). The error bars on the photometry have the same size of the symbols. The model fluxes 
and the COS spectra were convolved to have approximately a spectral resolution of 800, for display 
reasons. 
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(IPilbratd 



2010h . 



7. Conclusions 



On January 3rd and 5th 2010 we observed the planet hosting star WASP- 12 with the 
ESPaDOnS spectropolarimeter with the aim of looking for a global magnetic field, derive the 
stellar parameters and perform a precise abundance analysis. 

The WASP- 12 fundamental pa rameters an d iron abundance we obtained are in agreement 



with what previously published by 



Hebb et al. 



(l2009h : reff=6250+ 100 K, log ^=4.2+0.2, 



fmic=l-2+0.3kms~\ and [Fe/H] =0.32+0.1 2 dex. vs'mi is less than 4.6+0.5 km s~\ while fmacro 
lays within the range 4.75-7.0+0.6 km s~^ Dedicated observations to measure the amplitude of 
the RM effect are necessary to derive the real stellar v sin i. The resultant metallicity of WASP- 12 
is 2=0.021+0.002 dex. A detailed analysis of the HR diagram, with the use of isochrones and 
evolutionary tracks allowed to derive m ore accurate range s for the stellar age, mass and distance: 



the age of WASP-12 is between l.OGvr .Hebb et al 



and 1.49 M/Mq (the last value comes from 



Hebb et al 



(120091) an d 2.65 Gyr, the mass is between 1.23 



20091), and the distance to WASP-12 is 



between 295 and 465 pc. Our measurement of the Li abundance allowed us just to conclude that 
WASP-12 is older than 500 Myr. 

We performed a magnetic field search adopting the LSD technique revealing that the star 
does not show any magnetic field signature in Stokes V. A detailed analysis of the possible 
star-planet magnetic interaction w ould require a time-dependent an alysis of particular spectral 



lines, such as the Ca H & K lines (IShkolnik et al 



2003L 



2005L 



2008h and the Mg II UV resonance 



ines, but it has to be taken into account that WASP-12 shows a remarkably low stellar activity 



( Knutsonet al. 



20101) . that will be analysed in detail in a forthcoming paper. 



Given recent theoretical predictions (|Li et al 



2010|) and discoveries (IFossati et al 



20101), the 
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WASP- 12 system seems to be an ideal target to detect the presence of atmospheric pollution, 
due to the material lost by the planet. Therefore we looked for hints of pollution by looking for 
a correlation between the atmospheric element abundances and the condensation temperature, 
and by comparing the WASP- 12 abundance pattern with the one of other planet hosting and 
non-planet hosting stars, previously published by several other authors. Our analysis revealed just 
the presence of hints of atmospheric pollution, although only a differential analysis would allow 
to obtain firm evidences. One must also take into account the fact that it is not clear whether 



WASl P-12 would show the same kind of atmospheric pollution shown by the sun (|Melendez et al. 



20091) : the material coming from WASP- 12 b and falling onto the star is in a gas/plasma state and 
only a detailed modeling of the temperature and density structure of the accretion disk would 
show whether the material is condensing in dust grains. If dust grains are forming it is likely that a 
differential analysis would reveal the kind of pollution signatures obtained for the sun, otherwise 
all the material lost by the planet would fall onto the star, making then the pollution signature 
dependent to the unknown hydrogen content of the planet. 

With the use of the available HST calibrated spectra and of visible and infrared photometry, 
we looked for the presence of a circumstellar disk around WASP- 12, but without success. 
Probably high precision far infrared measurements, such as those possible with the Hershel 
satellite, may reveal the presence of a circumstellar disk. 
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